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We study the Lyapunov exponents of a two-dimensional, random Lorentz gas at low density. The 
positive Lyapunov exponent may be obtained either by a direct analysis of the dynamics, or by 
the use of kinetic theory methods. To leading orders in the density of scatterers it is of the form 
Aon Inn + Bon, where Ao and Bo are known constants and h is the number density of scatterers 
expressed in dimensionless units. In this paper, we find that through order (n^), the positive 
Lyapunov exponent is of the form ylonlnn + Bon + Aifi^lnn + Bifi^ . Explicit numerical values 
' of the new constants Ai and Bi are obtained by means of a systematic analysis. This takes into 

account, up to O(n^), the effects of all possible trajectories in two versions of the model; in one 
version overlapping scatterer configurations are allowed and in the other they are not. 
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I. INTRODUCTION 



The Lorentz gas is a model consisting of a set of scatterers that are fixed in space, together with a moving point 
, particle (or a cloud of mutually non-interacting point-particles) undergoing collisions with the scatterers. Here we 
will consider two variants of the two-dimensional version, where the scatterers are fixed hard disks, placed at random 
in the plane and the collisions of the point particle with the scatterers are elastic and specular. In the first version 
the positions of the scatterers are completely random, so different scatterers may overlap each other (this corresponds 
to the case of point scatterers with a moving particle of circular shape). In the second version the scatterers may 
not overlap each other, but each configuration satisfying this constraint has equal a priori weight (this corresponds 
to a hard-sphere interaction between the scatterers). The Lorentz Gas has proved to be very useful for studying 
, the general relations existing between dynamical systems theory and the non-equilibrium properties of many body 

■ systems. Explicit examples of such relationships encompass the escape rate formalism of Gaspard and Nicolis as 
I I well as the Gaussian thermostat formalism developed by Evans and co-workers 2, 3] and by Hoover and co-workers 

• The Lorentz gas is known to be chaotic due to the convex shape of the boundaries of the scatterers. Its 

' chaotic properties can be analyzed in fair detail, at least if the system is sufficiently dilute, in other words, the average 

■ distance between neighboring scatterers is large compared to their radius. For the model where the scatterers are 
•n ' placed on a periodic lattice (the Sinai billiard) Sinai has shown mixing and ei^odicity and demonstrated that 

, on large time and length scales the motion of the point particle is diffusive U Already much earlier Krylov 
(~| conjectured that to leading order the positive Lyapunov exponent is of the form nlnn, with h = na^ proportional to 
' ^ ] the density of scatterers. Subsequently Van Beijeren et al. 0, carried out kinetic theory calculations yielding 
(Ih ■ explicit expressions for the Lyapunov exponents of a disordered Lorentz gas (random configuration of scatterers) to 
' leading orders, i.e. up to 0{h) in the scatterer density. These studies confirmed Krylov's conjecture and provided 
explicit values for the constants a-nd i3o appearing in the leading terms A^fihifi -\- Bqu of the expansion of the 
largest Lyapunov exponent in terms of the scatterer density. The kinetic theory methods used for obtaining this 
result, employ an averaging over all allowed configurations of the scatterers. 

Up to quadratic order in the density of scatterers, the positive Lyapunov exponent can be calculated by similar 
kinetic theory methods, refining the averaging procedure described above, such that it (i) takes into account the effects 
of non-overlap of the scatterers (for the non-overlapping scatterer model) and (ii) accounts for the most important 
effects of correlated collisions. The additional contributions to the positive Lyapunov exponent resulting from (i) 
and (ii) yield the positive Lyapunov exponent correct through order n^. The purpose of this paper is to present 
a systematic analysis to calculate these additional contributions and to give explicit values for the new coefficients 
appearing in the expression Aonlnn -I- Bqu + Aiv? Inn -f BifP' for the positive Lyapunov exponent up to this order. 
This is done in two steps: in the first step, we calculate the positive Lyapunov exponent up to O(n^), assuming that 
the collisions between the point particle with the scatterers are uncorrelated. This embodies that the particle does not 
encounter any scatterer more than once, plus that all previous knowledge on the presence and absence of scatterers 
on the track of the point particle is ignored. However, in the non-overlapping scatterer model the influence of the 
non-overlap condition on the collision frequency is taken into account (for the overlapping model no corrections are 
needed for this at the given orders in n). In the second step, we consider (ii), and calculate corrections to the Lyapunov 
exponents due to correlated collision sequences in which the point particle either encounters the same scatterer more 
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than once, or its scattering probability is enhanced or suppressed by the knowledge resulting from previous collisions, 
or the absence thereof. This includes the immediate suppression of collisions as result of the non-overlap condition. 

The structure of this paper is as follows: in Sec. II, we briefly discuss the general theory of the Lyapunov exponents 
for a two-dimensional, random Lorentz gas. In Sec. Ill, we calculate the Lyapunov exponents through 0{n^) and 
obtain the constants Ai and Bi, defined above. We conclude the paper with some discussions in Sec. IV. We note 
here that the structure of the calculation presented in this paper is based on a thorough and intricate mathematical 
formalism, which is summarized in the Appendix. 



II. LYAPUNOV EXPONENTS OF THE RANDOM LORENTZ GAS IN TWO DIMENSIONS: GENERAL 

THEORY 



The random Lorentz gas consists of point particles of mass m, moving among a random array of fixed scatterers. In 
two dimensions, the scatterers are hard disks of radius a. We will consider two versions of the model: the Lorentz gas 
with overlapping scatterers, in which each configuration of scatterers a priori has equal probability (hence no weighting 
according to the amount of free volume left to the point particles!), and the Lorentz gas with non-overlapping scatterers, 
in which scatterers cannot overlap each other, but each configuration of scatterers satisfying this constraint a priori 
is equally likely. For a system with N scatterers in a two-dimensional volume V the number density of scatterers is 
n = N/V, and at low density n ^ 1. There is no interaction between point particles. These particles, therefore, travel 
freely during flights between collisions with the scatterers. The collisions between a point particle and a scatterer are 
instantaneous, specular and elastic. During flights, the equations of motion of a point particle are 

r = V = — , p = mil = 0. (1) 

m 

At a collision with a scatterer, the post-coUisional position and velocity, r_|_ and of the point particle are related 
to its pre-collisional position and velocity, r_ and , by 

■r+ = r_ , t)+ = — 2 ■ &) & . (2) 

Here & is the unit vector in the direction from the center of the scatterer to the point of collision (see Fig. 1). Notice 
that the dynamics generated by Eqs. (|1I2|I keep the speed of the point particle constant at v. The instantaneous 
velocity direction of the particle is given by v{t) = v{t)/v. 



FIG. 1: Collision between a point particle and a scatterer. 



The Lorentz gas with circular scatterers is a hyperbolic dynamical system. Due to the convex nature of the collisions 
with the scatterers, typically, the distance between two infinitesimally close trajectories in phase space, [r{t), v(t)] and 
[r + 6r{t),v{t) + dv{t)], increases exponentially with time. There are two non-zero Lyapunov exponents, which sum 
to zero. We denote the positive Lyapunov exponent by A+. Without any loss of generality, one can characterize the 
time evolution of the separation between two nearby trajectories by that of its projection onto v-space. The positive 
Lyapunov exponent then can be defined as 

A.= lim hm 1 m l^-fa, + f^)' . (3) 

The non-overlapping random Lorentz gas is generally supposed to be ergodic, even in the infinite-system limit (al- 
though we do not know of any proof of this). Therefore this definition of A_|_ should be (almost) independent of the 
choice of the initial point of the trajectory. For the overlapping Lorentz gas the situation is slightly more subtle: 
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in this model, if the volume becomes large enough, ergodicity will be broken. There will always be finite enclosures 
formed by three or more scatterers, from which a point particle cannot escape if it is trapped inside initially (see Fig. 
EJ. However, as long as the density of scatterers is below a critical percolation density, in the infinite system limit 
there will always be a single unbounded percolating region on which the motion of the point particles is diffusive on 
large time and length scales. We will be interested in the Lyapunov exponents of such particles only (even though 
particles moving inside an enclosure will also exhibit a positive Lyapunov exponent). 




FIG. 2: A situation for the Lorentz gas with overlapping scatterers, wherein the point particle cannot leave the enclosure 
formed by the four scatterers. 



As an alternative to Eq. Q, one can also choose to measure the separation of two nearby trajectories in r-space 
dSr{t) 

— after all, Sv{t) — — — — . In this representation the positive Lyapunov exponent becomes 

A,. „„ „„ il.,l^^|L±il. (4) 

T^oo |5r(to)HO T \5r(tQ)\ 

When calculating A+ using Eq. one may take advantage of the feature that the separation in velocity space 
between the two trajectories undergoes a change only at the collisions with the scatterers. If the point particle suffers 
k collisions between time and to + T, then [lO| 

A+ = hm hm AiVlnM^, (5) 

i—l 

where, Svi- and 6viJ^ are respectively the pre- and the post-coUisional separation in velocity space between the two 
trajectories at the i-th collision. 

On the other hand, to obtain the time evolution of |(5r(t)|, one may introduce another dynamical quantity, called 

the radius of curvature, and defined as p = v \^J'^P, ^7^^ — SIS El- It characterizes the divergence between 

\ov(t) — v[t) ■ dv(t)\ 

neighboring trajectories and in two dimensions it may be defined as the distance of the actual positions on a pair of 
such trajectories to their mutual intersection point (with p positive if this point is found in the past). In terms of the 
radius of curvature the positive Lyapunov exponent may be expressed as j5j , 

V /■*"+^ dt 

Am = hm — / -^T- (6) 

T^^ T J,^ Pit) ^ ' 

During a free flight, the equation of motion for p is given by Q 

p = v. (7) 

At a collision with a scatterer, the post-collisional radius of curvature /9+ is related to the pre-coUisional radius of 
curvature p- by 0,11,01 • 

p+ p- a cos 
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which is well-known from geometric optics. Here is the collision angle, i.e, coscf) = \v- ■ a\ = \v+ ■ <t| (see Fig. 1). 
Equations iQ and for a given initial condition [r{to), ■w(io)] and a fixed spatial arrangement of the scatterers, can 
be solved together to obtain p{t) as a function of r{to) and v{to)- Consider a set of trajectory pairs generated from the 
same reference trajectory by starting tangent trajectories with Sr{to) — and 6v{to) — Svq at a range of initial times 
to. One easily convinces oneself that the radius of curvature along the trajectory rapidly approaches a limiting value, 
p{r, v) = limt„^_oo pir, v, to), with p{r, v, to) the radius of curvature at the phase point (r, v) for a trajectory bundle 
starting out from [r{to),v{to)]- Combining Eqs. {T)) and one finds that as a function of decreasing to p{r,v,to) 
can be expressed as a rapidly converging continued fraction. In Eq. p{t) therefore can be replaced by p{r{t), v{t)), 
independent of initial conditions, and, assuming ergodicity, one may replace the long time average in Eq. © by an 
equilibrium average 0, S llO j . Besides an average over initial position and velocity of the point particle this involves 
an average over all allowed configurations of the scatterers. In the sequel we will denote it as [^.TsLllol 

. (J) . 

A very useful way of rewriting the positive Lyapunov exponent uses the relationship 

\Sv+\ _ p_ 



\6v^\ p+'' 



(10) 



with the p's defined again as the values on a trajectory starting in the infinite past. Combining this with Eq. ^ one 
finds that for an ergodic system the positive Lyapunov exponent may be expressed as 

A+=i.e/ln^\ , (11) 

\ P+f coll 

where fc is the average collision frequency, which is the inverse of the mean free time between collisions Tc, and the 
average now runs over the equilibrium distribution of collision configurations. Using Eq. (jSJ one may rewrite this 
further as 

A+=.c/ln "^"^^+,'^- \ . (12) 
\ «cos0 /^^u 



III. THE DENSITY EXPANSION OF THE POSITIVE LYAPUNOV EXPONENT 

We now present an analysis for calculating the positive Lyapunov exponent, A+, up to O(n^), inclusive. As 
mentioned before, this is done in two steps. In the first step, we assume that the collisions of the point particle with 
the scatterers are uncorrelated, implying that the probability density for a collision with collision angle 4> (see fig. 1) 
at all times is given as 

Pcoii (</<) = cos 0. (13) 

We calculate the resulting contribution to the positive Lyapunov exponent up to 0{n^), by means of both the velocity 
deviation method and the radius of curvature method and show that the two agree. 

In the second step, we calculate the corrections (SA+ resulting from the non-overlap condition (in the overlapping 
scatterer model) and from correlated collisions, i.e. either collisions of the point particle with a scatterer it has collided 
with before, since the initial time t = to, or collisions with a collision rate that is slightly enhanced by the available 
information on the absence of other scatterers on the preceding free path. 



A. Uncorrelated collision approximation 

In the velocity deviation method our starting point is Eq. 1121 First of all we note that the contribution 
j/c(ln(l/cos0))coii is given exactly for all densities by 

i'c(ln(l/cos^))coii = -i^c / — lncos(/) 

J-7T/2 2 

= i/c(ln2-l). (14) 
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So, it remains to calculate the contribution i^c (In '"'-'^ — . To investigate this, consider fig. El The collision 

\ a 

angles with the scatterers 1 and 2 are denoted as and respectively. From Eq.[7|it follows that the radii of curvature 
of 2 and 1 just before respectively just after the collision of the light particle are related as 

and from Eq. |S1 we find that 



pL'' =pi''+z;T, (15) 



a cos ( 

(1) ^ acos0 _ V 2 ^ 
^+ 2 d) a cos (/) ■ ^ 

^- 

^From the last equation we see that at low density of the scatterers p^"* is of order a and typically by an order of 




FIG. 3: A picture of a collision sequence. The solid and the dashed lines are two trajectories of the point particle, infinitesimally 
displaced from each other. Scatterer 1 and scatterer 2 are denoted by 1 and 2 respectively. 



n smaller than the free path vt. Therefore, in our calculations of Lyapunov exponents the contributions from 

in Eq. I|15|l only will contribute to higher density corrections and not to the leading terms. For similar reasons the 

last term in Eq. I|l()|) may be ignored completely in our calculations, as it will not contribute to corrections of order 

, 1 • /^ 2t;T + a(cos(/) + cosi/j) ' 

smaller than n . The term still to be calculated may then be written as Vci In ■ 



. This has 



coll 



to be averaged over the distribution of 0, -0 and r, which are all three assumed to be independent. The distribution 
of r, under the assumption of constant collision frequency, becomes a simple exponential. Therefore the resulting 
expression for the Lyapunov exponent becomes. 



A 



(U) 



7r/2 



-tt/2 

In^ 

ah'c 



, C0S( 
3 

2 



'^/^ , ,cosV 

^ JQ 



In 



2vT + a(cos 4> + cos ip) 



In- 



Av 



C 



(17) 



where C = 0.5772 • • • equals Euler's constant. To rewrite this in terms of the scatterer density, one has to express the 
collision frequency in terms of the latter. For non-overlapping scatterers the overall collision frequency depends on 
the configuration of scatterers, but on average it equals 



2nav. 



(18) 



for all densities. However, this is only true if one averages over all initial points for the point particle, that is, points 
inside the enclosures as well as on the percolating region. At low densities the collision frequency of a particle inside an 
enclosure will be much higher than the overall average value, so the collision frequency of a particle in the percolating 
area, which we are really interested in, has to be smaller. Fortunately, these corrections are at most of relative order 
compared to the leading term (because it takes at least three scatterers to make an enclosure), so in our present 
analysis they may be ignored. In the case of non-overlapping scatterers the collision frequency, for all densities, follows 
immediately from the ratio between the sum of the circumference of all the scatterers to the available free volume for 
the point particle, as 



(nov) _ 



2nav 



(19) 



6 



Substituting these expressions into Eq. (|17(l . we find that A'i'^-' for the respective cases behaves as 



^(U)(ov) 



A 



(U)(nov) 



2nav 
2nav 



1 

2fi 
1 

2fi 



1-C 
1-C 



2 

TT 



n 



1 



-7T, 3 In ■ 
2 V n 



In 2 
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3C 



(20) 
(21) 



In the radius of curvature method, is calculated by taking into account precisely the same dynamical aspects 
and approximations that have been used in the velocity deviation method above. The identity of the results may be 
established immediately by rewriting the integral in Eq. © as a sum of integrals between subsequent collisions of 
the point particle. By using Eqs. ijTJl and 115|) one immediately finds that the positive Lyapunov exponent may be 
expressed as 



X+ = Vci In ■ 



p+ . 



(22) 



coll 



which is obviously equivalent to Eq. Hll|l . 

The approximation that the collisions suffered by the point particle are uncorrelated is taken into account by the 
use of an extended Lorentz-Boltzmann equation (ELBE) — defined below — for the distribution function /(r, v, p, t) 
of the moving particle in (r, v, p)-space at low density of scatterers. This function describes the probability density of 
finding the moving particle at position r with velocity v and radius of curvature p, at time i, averaged over all allowed 
configurations of scatterers (remember that for given configuration of scatterers p is a uniquely defined function of r 
and v). Since the speed of the particle is constant we may replace v by the angular variable 9 describing the angle 
between v and the a;-axis. In equilibrium the distribution function f(v, 0, p) is a function of the radius of curvature 
only and the extended Lorentz-Boltzmann equation takes the form |a, [l3 



df{p) P cos(/) /acos0 A f°° , , J. f p'acos(f> A 

= acos^ + 2pJ /(/^ ) - '--/(Z^) ' (23) 

with 0(a;) the unit step function. In this case, the distribution function /(r, 9, p, t) satisfies the normalization condition 

/ dpf{r,9, p, t) = F(r, 9,t) = ^ (24) 
Jo 27r\/ 

Here F{r,9,t) is the distribution function of the point particle in {r,9) space. 

The quantity p' in the argument of the (5-function in Eq. H23|l , is the pre-coUisional radius of curvature that produces 
a post-collisional radius of curvature p. From Eq. H15|) and the fact that the free path vt for most inter-collision paths 
is of the order h^^, it follows that in most cases p' 3> a. As a result of this one may to leading order in the density 
simplify Eq. H23() by using the approximation [3l 1^ ll(| . 



p'acoscj) \ / acos( 



i0 + 2p' 



(25) 



The solution of Eq. (|23|l obtained with this simplification, will be denoted as /^"Hp)- The positive Lyapunov exponent 
to lowest order in the density of scatterers, A^"*, is obtained by using the distribution function /'■*''' (p) in Eq. ^ to 
calculate the ensemble average. We express the solution of Eq. f{p) as a power series expansion in h as 



fip) = f^'Hp) +Sf^°\p) + ■■■ (26) 



and obtain the expression of A^^ as 



(27) 



We specifically point out here that to calculate (5/'*''(p), one cannot use the approximation in Eq. (|25|l any longer. 
The function Sf''°^p) may be calculated 0| by applying a successive approximation scheme to Eq. H23|l . Density 
corrections to Aq, contained in A^\ result from both f^^\p) and 5f^'^\p). This procedure for calculating these 



7 



corrections is algebraically more cumbersome than the velocity deviation method. Nevertheless, it equally leads to 
the results of Eqs. H2UI21|I . However, for systems that are not, on average, spatially homogeneous (e.g. a Lorentz gas 
with open boundaries) the velocity deviation method becomes more complicated and in fact the ELBE method is to 
be preferred [T^. 

In the second step of our calculations we consider corrections due to correlations between collisions. From Eqs. 
H12|l . (|15|l and H16|l it becomes clear that there are two types of corrections that have to be accounted for: 

• The approximation of /7_ by vt + a cos 0/2. 

• Deviations from the exponential Vc exp(— j/ct) of the distribution of the free flight time r between collisions. 

Through order the first type of corrections are only important in case vt is comparable to a, that is if two scatterers 
with which the light particle collides are close to each other. If this is the case, the light particle may recoUide an 
arbitrary number of times with cither scatterer and at each of the collisions the actual value of p_ has to be used 
instead of the uncorrelated collision approximation. This will be discussed further in the next subsection. The second 
type of corrections to order tt? only occur for the case of non-overlapping scatterers. They will be investigated in the 
next-to-next subsection. 

B. Corrections due to nearby scatterers 

In case two scatterers have a mutual distance comparable to their radius a, a first collision of the point particle 
with either of them with fairly large probability will lead to recollision sequences such as the one shown in fig. ^ 




FIG. 4: An example scattering sequence involving multiple collisions with two scatterers. 

For the first two collisions, at P and Q, the approximation of p- by ur + o cos with r and (p the preceding free 
time respectively the collision angle at the previous collision, is still adequate. For all subsequent collisions p- has to 
be calculated from Eqs. ^ and ((T^ . This gives rise to a correction 

) = .7ln , 2p +acos^ \ 
\ 2vT + a(cos ip -h cos 0) / 

to the positive Lyapunov exponent. The average has to be taken over all two-scatterer configurations with nearby 
scatterers and over all points on the circumference of scatterer 1 where a collision may occur (so in the overlapping case 
points covered by 2 are excluded) and over all allowed values of the coUison angle -0, weighted again by cos ■0/2. In fact 
only those parameters that correspond to a recoUison with 1 give rise to non-vanishing contributions at the relevant 
orders of density. For the overlapping scatterer model the initial configurations of the scatterers do include overlapping 
ones; for the non-overlapping model such configurations are excluded. We evaluated the averages numerically, with 
the results 

^^(ris)(ov) ^ _o.38(2naw)n and 

^^(ns)(nov) ^ -Q,2S{2nav)h. (29) 

Notice that the probabilities that these collision sequences are interrupted by other scatterers can clearly be ignored 
at the orders of the scatterer density we are considering presently. 
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FIG. 5: For given scattering angles (j) and %p the free path length between subsequent collisions has a minimal value b{<j),tp). 
This minimum is determined by the condition that the two scatterers touch each other. 




FIG. 6: After a long free flight the probability that a collision will occur like the one shown with scatterer 2, is slightly enhanced 
by the fact that no scatterers like 1 can be in the way. If the collision angles are ip and (j) respectively, the distance over which 
the absence of 1 influences the probability of finding 2 is just b{(l),tp). 



C. Corrections to the free-path length distribution 



The assumption of a constant collision frequency, independent of the past, is not entirely correct. For the case 
of overlapping scatterers this only becomes manifest when one considers corrections of relative order ri^ . So for our 
present purposes we may ignore this. For non-overlapping scatterers corrections already show up at the order h. 
These appear in two forms: shortly after a collision the probability for a next collision is reduced by a shadowing 
effect: many scatterer locations that would give rise to such a collision are forbidden by the non-overlap condition. 
This is compensated by an anti- shadowing effect at long times; the absence of scatterers in a strip of width 2a around 
the free path of the point particle increases the probability of finding a scatterer with which this particle will collide. 
Obviously, the combined influence of these two effects on the overall collision frequency has to vanish, but there is 
a shift of collision probability towards longer times, shifting the free-path length distribution equally to somewhat 
larger values of the free path. 

The shadowing effect is illustrated in fig. |S| For subsequent scattering angles (j) ^-nd ip the minimal path length 
between the collisions is given by 



&((/>,■)/')= a ■\/4 — (sin cj) — sin -0)2 — cos (f> — cos ip . (30) 

The contributions from shorter path lengths have been counted erroneously in and should be subtracted. 

This gives rise to a correction of the form 

af^'^)^-' = -{2navf j\d4>'^ |"^dV.^yJ|''''''"'rfrln(^+cos<^ + cosV^), 

= -0.0976 (2nat;)n. (31) 



The anti-shadowing effect occurring for long free times, is illustrated in fig. For a collision with collision angle 
there is with certainty no obstruction from a scatterer that would have given rise to a previous collision with scattering 
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angle (f) within a preceding distance 6(0, ^Z"). This enhances the colhsion frequency by an amount of 

7r + 18V3- 32 



, cos ( 



dtp — --^ nab{(f>,'ilj) 



24 

0.0966 • • • i^c 



(32) 



One should also take into account, however, that the probability for a collision at long times is enhanced somewhat 
by the initial suppression of the collision frequency through the shadowing effect. The combined effects lead to a long 
term collision probability density of the form 



p,,,(t) = (^, + 5^0t)),-K+*.<"')(*^*.). 
Here td is the average delay time due to the shadowing effect. It is given by 



td 



1 



and related to Si^c^^ by 

The ensuing correction to the positive Lyapunov exponent now may be obtained as[l 

> (as) (nov) 



(33) 



(34) 



(35) 



In — -C 
2n 



Collecting the contributions to from Eqs. (|20|l . (|21|l . (|29|) . H31|) and H36|l . we obtain our final result. 



S\ 



(ov) _ 



dX 



(nov) 

+ 



2nav 



2nav 



In — + l- C + -n ln- C - 0.242 



2h 2 \ h 6 



In h 1 - C + -n 31n In 2 

2h 2 \ fi 



and 



II 



3C - 0.240 



I8\/3-32 /, I 



I27r 



In C 

2n 



(36) 

(37) 
(38) 



IV. DISCUSSION 



In this paper, we have calculated the Lyapunov exponents of a two-dimensional, random Lorentz gas at low densities 
up to Oiji?) in the density of scatterers. This calculation was carried out in two parts: in the first part we assumed 
that subsequent collisions between the light particle and one of the scatterers are uncorrelated. In the second part we 
calculated the effects of correlations between collisions and, in the case of non-overlapping scatterers, those resulting 
from the non-overlap condition. The effects of repeated recoUisions of the light particle with two scatterers that are 



close to each other, were calculated numerically, as well as the contribution from 5\ 



(shad) (nov) 
-I- 



AU other contributions 



to the positive Lyapunov exponent were obtained analytically. For the sake of brevity and understanding, we have 
presented the method in this paper in fairly intuitive terms, as opposed to the more formal mathematical structure 
upon which all calculations have been based originally A short summary of this formalism, is included in the 
Appendix at the end of this paper. 

We argued that our analysis covers all density correction to the Lyapunov exponent through order . We have no 
rigorous proof of this, but it is easy to support this claim by simple power counting arguments. 

The formalism was developed in the context of calculating the Lyapunov exponents of a two-dimensional Lorentz 
gas at low densities up to 0{p?), but the method itself is not limited to the calculation of the Lyapunov exponents 
alone - it can be used for a density expansion of other dynamical quantities of a two-dimensional Lorentz gas as well. 
We note that the formalism as well as the somewhat intuitive analysis in this paper, in principle can be extended 
so as to obtain expressions for the Lyapunov exponents up to higher orders in the density of scatterers. Though 
it would be very desirable to have reliable results for all densities, we expect a systematic extension of the present 
methods will be very laborious and complicated. Perhaps it will be possible to find non-systematic, but still accurate 
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approximations, such as Enskog equations or ring kinetic equations in the kinetic theory of transport processes. In 
that case the present theory may be very valuable in providing guidelines and testing criteria at moderate densities. 

A more modest goal would be an extension of the present analysis to systems out of equilibrium. Here one may 
think first of all of systems with escape and systems under the action of a driving field combined with a gaussian 
thermostat. This should certainly be doable. 
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APPENDIX 

Here we want to illustrate how the results obtained in this paper may be obtained in a systematic way from a 
diagrammatic analysis of the dynamics of the Lorcntz gas. 

Startin g po int is the binary collision expansion (BCE), explained for the Lorentz gas, respectively the hard sphere 
gas in [TtiIISII . The dynamics of the point particle is expressed as a series of convolution products involving sequences of 
alternating free streaming and binary collision operators. The free streaming operators describe the ballistic motion of 
the point particle between two collisions with a scatterer. The binary collision operators consist of a real and a virtual 
part. Both of them check for the conditions of a collision to be satisfied. The real collision operator implements the 
velocity change of the point particle, as given by Eq. JSJ. The virtual collision operator leaves the velocity unchanged 
but multiplies by —1. It always produces corrections to simpler terms in the binary collision expansion, which ignore 
the collision and let the point particle move freely through the scatterer. 

The binary collision expansion is built by starting from a free streaming operator ignoring all collisions occurring 
in reality. The first corrections consist of terms containing one collision operator, with the real one generating a 
path that contains one collision with a specific scatterer and the virtual one subtracting the free streaming terms for 
all initial configurations from which free streaming indeed gives rise to a collision with that given scatterer. Next, 
two-collision terms give rise to paths with two real collisions and correction terms with either one real and one virtual, 
or two virtual collision terms, and so on. 

The non-overlap condition between point particle and scatterers may be accounted for at the initial time, but, as 
shown by Van Beijeren and Ernst it has great advantages to postpone this, as much as possible, until the times 
of the first collision of the point particle with any specific scatterer. This may be done efficiently by choosing, in each 
term of the BCE, the first collision operator with any given scatterer, unless it is the last collision with this scatterer 
as well, to be a T operator, rather than a T operator, used for all remaining collisions. The difference between these 
two binary collision operators is that the T operator counts virtual collisions at the instant that the point particle 
leaves the scatterer, whereas the T operator counts them at entrance. 




FIG. 7: Cancellations between real and virtual collision diagrams in binary collision expansion for repeated collisions between 
the point particle and two scatterers separated by distances of 0{a). 

The uncorrelated collision approximation is reproduced by keeping in the BCE exactly those terms in which each 
scatterer appears at most once. A good way of rearranging these terms is separating each of the binary collision 
operators into its real and virtual part, and taking together all events that have the same sequence of real collision 
operators. Terms with arbitrary number of virtual collisions between two subsequent given real collisions occurring 
at times tj and ij+i, sum to yield the damping term exp[— fB(ij"-i-i — tj)], with = 2nav the low density limit of the 
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collision frequency. For the case of non-overlapping scatterers one may extend the analysis so as to replace the low 
density collision frequency vb by the full collision frequency Vc- To this end one has to replace the collision operators 
in the BCE by sums of operators representing the product of the collision operator with the pair correlation function 
of the point particle and a scatterer, at contact. The value of the latter just equals (1 — 7rn)~^, thus reproducing Eq. 
(|19|1 . Notice that also the particles responsible for the static correlations at collisions occur in these BCE terms at 
just one single point in time. 

Next, consider correlations due to the same scatterer being present at more than one collision. As in the main text, 
we distinguish between nearby cases, where all scatterers involved are within mutual distances of the order a, and 
events in which at least one of these distances is of the order of the mean free path, £ = a/ {2n). In the nearby case each 
additional scatterer reduces the weight by an additional factor proportional to n, so in the approximation considered 
here, we may restrict ourselves indeed to events involving just two scatterers. As stated in the main text, both for 
the overlapping and for the non-overlapping case we obtain contributions from events in which there is at least one 
recollision with either of the two scatterers. As an example, consider an event of the same type as exposed in Fig. 21 
but now with the collisions preceding virtual collisions at P or Q added (see fig.E)). Consider first the contributions 
to Eq. (|l()(l from collisions at R. For all BCE-events containing a collision at S one has exact cancellation from the 
contribution with a real and that with a virtual collision at S. So we only need consider events for which the collision 
at R is the last one. The event with just one preceding collision, at Q, was included already at the uncorrelated 
collision approximation. The event where a virtual collision at P is followed by real ones at Q and R subtracts the 
contribution from the event with subsequent real collisions at T, Q and R in the uncorrelated collision approximation. 
The contribution with subsequent real collisions at P, Q and R finally, adds the actual contribution from the collision 
at R. Next, consider the contributions from the collision at S. The event with a virtual collision at Q subtracts 
the uncorrelated event with collisions at U, R and S. The events with a virtual collision, respectively no collision 
at P, followed by real collisions at Q, R and S, cancel each other and finally again, the event with all collisions real 
gives the actual contribution. These reasonings are easily generalized to collision sequences of arbitrary length. They 
immediately lead to the result of Eq. 




FIG. 8: Limitation of the collision angle to 0{a/£) for recollision events involving multiple scatterers that are separated by 
distances of 0{£). Notice that the trajectories between scatterers 1 and 3 are confined within the shaded grey area, which limits 
the variation of the collision angle at R within a range A(f) ^ a/£. 



In the non-overlapping case the contributions from consecutive collisions with mutually overlapping scatterers 
need special consideration (in the overlapping case such configurations play no special role). In the diagrammatic 
representation of the BCE the overlap gives rise to a non-vanishing Mayer-bond between the overlapping scatterers 
and collision events like the one shown in Fig. |31 correspond to sections of diagrams with two different T-operators 
connected by a Mayer-bond between the two scatterers involved. Either of the two collisions may be real or virtual. 
The case of two real operators leads to the correction ^^(^had)(nov) subsection IIII Bl The case of a real collision 
followed by a virtual one is responsible for the delay time td defined in Eq. H34|l . A virtual collision followed by a real 
one is responsible for the enhancement Svc^"^ of the collision frequency as specified in Eq. H33|) . The diagram segment 
with two virtual collisions accounts for the enhancement by Sv^^^ of the damping factor in the survival probability 
during free fiight. Finally, events involving recollisions between overlapping scatterers cancel exactly against the same 
events without the Mayer-bond between the two scatterers included. In subsection IIII Bl this was accounted for by 
restricting the contributions from recollision events involving two nearby scatterers to non-overlapping configurations 
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in the non-overlapping model. 




FIG. 9: Cancellations between real and virtual collision diagrams in binary collision expansion for repeated collisions between 
the point particle and multiple scatterers separated by distances of 0{£). 

Recollisions where the point particle returns to a scatterer from a distance of the order of the mean free path, 
restrict the allowed collision angle at the preceding collision to an angular range of the order a/£ ^ n (see Fig. 
Therefore, at the level of corrections to A_|_ restricted to order n^, no further restrictions may be imposed on the free 
flight lengths and collision angles of the intermediate collisions and all of these lengths typically are of order £. Now 
consider, as an example, contributions to (5A+ resulting from the recoUision at S in Fig. According to Eq. 111(1 

these come from the averages of In — at S, in case the recoUision is real, or at T, in case it is virtual. In either 

P+ . . . . 

case, the contribution in which the collision at P is virtual, cancels that where the latter is real, through the order 
of n considered. The same holds true for the contributions from all later collisions. Collisions before the recoUision 
of course are not influenced by it at all. Hence, we conclude that corrections to A+ from recoUision events involving 
distances of the order £ are absent through order and we were justified in ignoring these in subsection IIII Cl 
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